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Abstract 

After  deriving  a  set  of  dynamic  equations  governing  the  dynamics  of  a  Tethered  Satel¬ 
lite  System  (TSS),  stabilizing  tension  control  laws  in  feedback  form  are  derived.  The  tether 
is  assumed  rigid  and  massless,  and  the  equations  of  motion  are  derived  using  the  system 
Lagrangian.  It  is  observed  that,  to  stabilize  the  system,  tools  from  stability  analysis  of 
critical  nonlinear  systems  must  be  applied.  This  paper  employs  tools  related  to  the  Hopf 
Bifurcation  Theorem  in  the  construction  of  the  stabilizing  control  laws,  which  may  be 
taken  purely  linear.  Simulations  illustrate  the  nature  of  the  conclusions,  and  show  that 
nonlinear  terms  in  the  feedback  can  be  used  to  significantly  improve  the  transient  response. 


I.  Introduction 


The  topic  of  Tethered  Satellite  Systems  (TSS)  has  received  considerable  attention 
in  recent  years  (e.g.,  [l]-[8],  [13],  [14]).  Potential  applications  of  these  systems  include 
deployment  and  retrieval  of  satellites,  aiding  in  space- assembly  tasks,  use  of  electrodynamic 
tethers  for  electric  power  generation  [12,  p.  3-29],  and  tethering  platforms  above  the  Space 
Station  for  observing  stellar  and  planetary  objects  [12,  p.  3-85].  For  other  potential 
applications,  the  reader  is  referred  to  Rupp  and  Laue  [1]  and  the  NASA  publication  [12]. 

In  this  paper,  we  focus  on  the  issue  of  stabilization  of  a  tethered  satellite  system  during 
the  station- keeping  mode.  Specifically,  consider  the  TSS  depicted  in  Figure  1.  Here,  a  large 
satellite  is  tethered  to  a  smaller  subsatellite,  and  the  configuration  is  in  a  circular  orbit 
around  the  Earth.  During  station-keeping,  a  subsatellite  is  controlled  so  as  to  follow  a 
prescribed  orbit  to  within  a  set  tolerance  [15,  p.  220].  By  assuming  the  subsatellite  to 
be  much  more  massive  than  the  satellite,  and  that  the  satellite  follows  a  perfect  circular 
orbit,  we  are  able  to  focus  attention  on  the  station- keeping  control  of  the  subsatellite. 
This  is  accomplished  through  the  design  of  tether  tension  control  laws  in  feedback  form 
which  result  in  regulating  the  position  of  the  subsatellite  relative  to  the  satellite,  while 
simultaneously  regulating  the  tether  length  at  a  prescribed  nominal  value.  The  proposed 
tension  control  law  is  implemented,  say,  using  a  reel-type  mechanism. 

The  paper  makes  use  of  several  simplifying  assumptions.  For  instance,  the  tether  is 
assumed  rigid  and  massless.  With  these  assumptions,  the  TSS  can  be  described  by  a  set  of 
ordinary  differential  equations.  The  system  Lagrangian  is  used  to  obtain  these  equations. 
Next,  we  observe  that  linear  feedback- type  tension  control  laws  can  place  all  but  two  poles 
of  the  system.  These  two  poles  are  a  complex  conjugate  pair  of  pure  imaginary  eigenvalues 
of  the  system  linearization.  To  stabilize  the  system,  therefore,  tools  from  stability  analysis 
of  critical  nonlinear  systems  must  be  applied.  Our  approach  is  to  use  the  technique  of 
[10],  in  which  Hopf  bifurcation  calculations  are  employed  to  construct  stabilizing  control 
laws.  First,  a  class  of  purely  linear  stabilizing  feedback  control  laws  are  given.  Next, 
nonlinear  stabilizing  control  laws  are  developed.  Simulations  are  presented  which  allow 
one  to  compare  the  transient  response  of  the  system  with  the  two  types  of  feedback.  The 
additional  freedom  afforded  by  the  inclusion  of  nonlinear  terms  can  be  used  to  obtain  a 
significant  improvement  in  the  speed  of  the  transient  response. 


2 


Notation 

E  -  Earth 
S  -  Satellite 

m  -  Subsatellite  and  subsatellite  mass 
G  -  Gravitational  constant 

M,  ms  -  Mass  of  the  Earth,  mass  of  the  satellite 

(®m ,  Dm i  Zrn)  -  Earth-based  rotating  Cartesian  coordinates  of  subsatellite,  with  zm  in  the 
local  outgoing  vertical  direction,  and  xm  in  the  direction  of  motion  of  the  subsatellite 
in  its  orbit  (see  Figure  2) 

(x a,ys,zs)  -  Earth-based  rotating  coordinates  of  the  satellite 
-  Inertial  coordinates  of  subsatellite 
(&s,ys,  %s)  -  Inertial  coordinates  of  the  satellite 
0,  -  Constant  angular  velocity  of  the  satellite  in  circular  orbit 

6,  (j>  -  In-plane  angle  and  out-of-plane  angle  of  subsatellite  relative  to  local  vertical 
uitf,  <j),  us  Q,  t  -  Tether  length,  v  £ 

r0,  rm  -  Radius  of  the  satellite  orbit,  radius  of  subsatellite  orbit 
TO-,  T<t>  -  Generalized  torques  in  directions  6 ,  (f> 

Ft  -  Generalized  force  along  tether 

F  :=  e  cos  cj>^  +  +  Fel,  where  a  hat  indicates  a  unit  vector  in  the  given  direction 

Re(-)  :=  Real  part  of  (•) 

Jto(-)  :=  Imaginary  part  of  (•) 

i  := 

II.  System  Model 

Referring  to  the  depictions  in  Figures  1  and  2,  a  mathematical  model  of  the  TSS  may 
be  derived.  Assume  the  satellite  and  the  subsatellite  are  point  masses  and  the  tether  is 
massless  and  rigid.  Moreover,  take  the  mass  of  the  satellite  to  be  much  larger  than  that 
of  the  subsatellite  (i.e.,  ms  m).  This  allows  us  to  take  the  center  of  mass  of  the  TSS 
to  be  the  satellite,  and  to  consider  the  satellite  as  being  in  a  circular  orbit  around  the 
Earth.  In  addition,  the  gravitational  attraction  between  the  subsatellite  and  the  satellite 
is  neglected. 

It  is  evident  from  Figure  2  that  we  have  the  relationships 

xm  =  £  cos  <j>  sin  9  (1) 
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Vm  —  £  sin  <j) 

Zm  —  J’O  +  t  cos  <f>  COS  9 

rm  —  r  o  +  +  2roi  cos  <f>  cos  9. 


Also, 


cos  tit  0  sin  tit 
0  1  0 
-  sin  tit  0  cos  tit 


xs\  (  cos  tit  0  sin  tit  \  /  0 

Vs  =  o  1  0  0 

zs  )  \  —  sin  tit  0  cos  tit )  Vro 


where  the  equations  above  fix  a  particular  choice  of  time  reference. 

Since  the  tether  is  assumed  to  be  massless,  the  total  kinetic  energy  of  the  system  is 
given  by 

1  ,  x  2  x2  %2.  1  -x2  x2  x2  x 

KE  =  —ms(xs  +y9  +  zs)+  -m(xm  +  ym  +  zm) 

=  i msti2r1  +  ~m{£2  +  £24>2  +  £?  cos2  </>(#  +  ft)2  +  ti2r2 
+  2tir0icos(f>sm9  —  2tiro£  sin  <j)  sin  9<j>  +  2tir0£  cos  (j>  cos  9(9  +  O)}.  (5) 


The  potential  energy  of  the  system  arises  solely  from  gravity  and  is  given  by 


PE  = 


GMrris  GMm 


Moreover,  since  the  satellite  is  assumed  to  be  in  a  circular  orbit,  it  is  in  a  zero-g  orbit. 
Thus 


GMms 


=  msti2r0, 


or  more  succinctly  GM  —  ti2r%.  Writing  the  expression  for  the  system  Lagrangian  L  = 
KE  —  PE  and  invoking  the  Lagrangian  formulation  of  dynamics,  the  dynamic  equations 
of  the  system  are  found  to  be 


p  .  . 

tq  =  m£2  cos2  4>{9  +  2 -($  +  ti)  —  2tan<^>(#  +  ti)(j> 
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(6) 


n2r0  sin#  _  r|_ 

'  *  /V1  3  /J 


£  COS  (j) 


£  : 


r,£  =  +  2-<^>  +  cos<^sin<^(0  +  Q)2 

+  ^p-  cos  ^  sin  <^(1  -  ^f)} 


(7) 


F>  =  m{t?  —  ^(^)2  —  £  cos2  ^(0  +  fi) 


m 

2  j/d  i  o\2 


+ 


f}2r^ 


02ro  cos  (j>  cos  9(\ - — )}. 


(8) 


For  the  limiting  case  r0  £,  we  have  rm  ~  r0 ,  and  Eq.  (4)  implies 

-3  £ 


1 - ~  3  cos  (j>  cos  0- 


ro 


(9) 


Using  Eq.  (9),  the  approximate  motion  equation  of  the  system  for  the  case  ro  >  t  is  found 
to  be 

F  =  m£{£  —  £&' 2  —  £  cos2  <j)(9  +  f))2  +  £Q?  —  3 £l2£  cos2  4>  cos2  9 } 

+  m9{9£  cos  <f>  +  2(6  +  Q)(£  cos  (f>  —  £<j>  sin  <j>)  4  Z£9?  cos  9  cos  (f>  sin  9} 

+  m(f>{£(f)  +  2£<j)  4-  £  cos  (j)  sin  <f>($  +  fi)2  4-  3^fl2  cos2  9  cos  <f>  sin  (j)}, 

which  agrees  with  the  model  derived  by  Arnold  [2]  using  the  gravity  gradient  method. 
Note  that,  in  the  analysis  of  the  following  sections,  we  do  not  assume  r0  >  £. 

III.  Analysis  and  Control  in  the  Station-Keeping  Mode 

III.l.  Model  in  State-Space  Form 

Suppose  cos  <j>  ^  0  (i.e.,  </>  ^  ±f )  and  let  the  applied  tension  force  be  the  only  external 
force  acting  on  the  system.  Thus,  for  instance,  we  neglect  effects  of  a  rotating  atmosphere, 
the  Earth’s  magnetic  force,  solar  radiation,  and  the  oblateness  of  the  Earth.  We  do  not 
take  into  account  the  mechanism  for  commanding  the  desired  tether  tension,  although  one 
can  imagine  it  to  be  controlled  by  a  reel  mechanism. 

Since  (i)  we  have  modeled  the  satellite  and  subsatellite  as  point  masses,  (ii)  the  tether 
is  assumed  rigid,  and  (iii)  there  are  no  external  forces  besides  the  commanded  tether 
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tension,  we  conclude  that  the  generalized  forces  acting  on  the  subsatellite  are  Ft  —  T, 
Tg  =  0,  =  0.  Eqs.  (6)-(8)  can  now  be  rewritten  in  state  space  form  as  follows: 

i  =  (10) 

2v  1  Q2rn  r3 

sin(2^)(u>^  +  ft)2 - —  cos  9  sin  <j>{\ - |-)  (11) 

0=U6  (12) 

2v  .  ft2rosin0.  .  .  . 

^  — (u>0  +  ft)  +  2tan<£(u;0  +  ft)u^ - - - — (1 - — )  (13) 

-t  b  COS  (p  ^ m 


i  =  V 

V  =  +  P  cos2  <£(W0  +  ft)2 - 

?  m 

r3  T 

+  ft2ro  cos  6  cos  ^(1 - 1~)  H - , 

rt  m 


(14) 


(15) 


For  the  case  in  which  the  tether  length  is  held  constant  (i.e.,  P  =  v  —  0,  P  =  P*  =  a 
constant),  the  conditions  for  an  equilibrium  point  are  6  —  nn  and  (f>  =  rmr,  where  n,  m 
are  integers.  (We  disregard  another  apparent  possibility  for  <j>,  since  the  corresponding 
equation  has  no  solution  for  (j>  when  0  =  mr.)  At  the  equilibrium  point  (0, 0, 0, 0)  when 
only  Eqs.  (10)-(13)  are  considered  with  v  =  0,  the  linearized  system  of  Eqs.  (10)-(13)  has 
the  two  conjugate  pairs  of  pure  imaginary  eigenvalues 


where 

rm,o '■=  r0  +  P* ■  (18) 

The  eigenvalues  Ai,2  are  associated  with  Eqs.  (10)  and  (11),  i.e.,  with  the  out-of-plane 
dynamics,  and  A3j4  are  associated  with  Eqs.  (12)  and  (13),  i.e.,  with  the  in-plane  dynamics. 

The  appearance  of  pure  imaginary  eigenvalues  suggests  the  possibility  of  oscillations 
near  the  equilibrium  point  (0,0, 0,0).  Specifically,  if  the  reel  mechanism  acts  like  a  latch, 
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resulting  in  a  fixed  tether  length,  the  system  may  have  librations  with  two  distinct  fre¬ 
quencies  along  with  the  orbital  motion. 

In  the  sequel,  we  consider  the  problem  of  designing  tension  control  laws  rendering  the 
TSS  asymptotically  stable  in  the  station-keeping  mode.  The  main  difficulty  will  be  the 
presence  of  the  two  pairs  of  pure  imaginary  eigenvalues,  and  the  uncontrollability  of  one 
of  these  pairs. 

The  conditions  for  an  equilibrium  point  of  system  (10)-(15)  are 


0  =  02  sin  4>{  cos  (f>  +  -y  cos  0(1 - ~)  } 

^  rm 

n  _  sinfl  _  r|_x 
i  cos  4> 

O  2r3f  r3  rp 

0  =  £Q,2  cos2  <j>  -  ■ — +  tl2r0  cos  8  cos  (j>(  1 - ~)  H - 

rl  rt  m 


(19) 

(20) 
(21) 


where  the  applied  tension  force  T  (applied  through  a  reel  mechanism)  may  be  constrained 
to  satisfy  (21).  From  the  definition  of  out-of-plane  angle  <j>,  we  have  -f  <  <j>  <  f .  There 
are  hence  only  two  equilibrium  points:  (0,  0,  0,  0,  £*,  0)  and  (0,  0,  ir,  0,  £* ,  0)  if  the  tether 
length  is  fixed  at,  say,  t  =  t* .  In  our  paper  [9],  it  is  observed  that  the  set  <j>  =  0,  =  0 

is  an  invariant  manifold  for  Eqs.  (10)-(15),  regardless  of  the  form  of  the  tension  control 
law  T.  Although  this  implies  that  the  system  (10)-(15)  is  uncontrollable,  we  find  below 
that  there  does  exist  a  control  strategy  stabilizing  the  system.  With  an  assumed  rigid  and 
massless  tether,  there  may  appear  to  be  no  constraint  on  the  value  of  the  applied  tension 
force  T.  In  reality,  however,  the  tether  is  not  rigid.  Thus  the  subsatellite  cannot  be  pushed 
away  from  the  satellite  by  the  applied  tension  force  through  the  tether.  Note  that,  in  this 
paper,  the  sign  convention  implies  that  a  positive  value  of  tension  would  correspond  to  a 
slack  tether.  Hence,  in  the  sequel  we  restrict  the  applied  tension  force  T  to  assume  only 
nonpositive  values.  Although  the  conclusions  we  will  reach  will  also  apply  to  the  model 
(10)-(15)  in  the  absence  of  this  restriction,  they  would  no  longer  relate  to  the  physical 
problem. 
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III. 2.  Stabilization  for  In-Plane  Angle  Near  0  =  0 

Let  Xo  =  (0,0,0,0,£*,0)T  and  X  =  x  —  x0,  where  x  =  (<f>,  u 6,  u>$,  £,  v)T.  Then  the 
Taylor  expansion  of  the  right  side  of  Eqs.  (10)-(15)  is,  to  third  order  in  X, 

4-X  =  LoX  +  Qo(X,X)  +  Co(X,X,X)  +  eU  +  e—  (22) 


where  the  matrix  Lo,  the  quadratic  form  (Jo,  the  cubic  form  Co,  the  vector  e  and  the  scalar 
U  are  given,  in  terms  of  parameters  a\  defined  in  Appendix  A,  by 

/  0  1  0  0  0  0\ 

-a?  0  0  0  0  0 

r  _  0  0  0  1  0  0 

L°  ~  0  0  -a\  0  0  -§?• 

0  0  0  0  0  1 

\  0  0  0  2£*Ll  a3  0  / 

0 

-2 £l(j>LO0  +  a.12 <f>£  - 

0 

ai29£-  j?loqv  +  2Q,<j>o;fj>  +  j^£v 

0 

04  02  +  £* <jOq  +  2QuJg£  +  0,5  4>2  T  £* T  0\3£2 

0 

a%92§  —  uj2e(j)  +  a7(j)3  +  a8cj)£2  + 

0 

a903  +  ai09(j)2  +  a%9£2  +  2<j)u^,ujg  +  j^uglv  -  -Pv 

0 

o,i4.92£  -|-  w“q£  —  2 £*£lu)Q(f)2  a,\\<j)2£  -\-  T  Qis^3 

e  =  (0,0,0,0,0,1)T 

_  (3rgr  +  3r0l*2  +  l*3)ft2 

(r0  +  £*)2 

where  £■.=  £  —  £*.  The  expressions  above  have  been  verified  using  the  code  MACSYMA.1 
Case  1:  Linear  State  Feedback 

1  MACSYMA  is  a  trademark  of  Symbolics,  Inc.,  Cambridge,  MA. 
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Our  first  design  is  that  of  a  tension  control  law  in  linear  state  feedback  form  which 
stabilizes  the  system.  The  design  is  carried  out  in  two  steps.  The  first,  addressed  in  Lemma 
1  and  Corollary  1  below,  is  to  give  conditions  on  the  linear  state  feedback  ensuring  that 
four  of  the  eigenvalues  of  system  (10)-(15)  are  moved  to  the  left  half  of  the  complex  plane. 
These  eigenvalues  correspond  to  the  “in-plane  variables”  (6,  ujq,  £,  v). 


Lemma  1.  If  the  following  conditions  hold,  then  the  tension  control  force  T  =  m(—U  — 
k\0  —  k2u>e  —  k3£  —  k±v)  stabilizes  the  “in-plane  Jacobian  matrix”  at  the  equilibrium  point 
xo ,  i.e.,  the  Jacobian  matrix  of  Eqs.  (12)-(15)  with  respect  to  the  vector  (<?,u>0,£,  u): 

(i)  ki  >  0 

(ii)  bub2,b3  >  0 

(iii)  &4&1&2  —  b\  —  k\b3  >0 


where 


h  =  k. 


(2  r0£*2  +  £*3)Q2  2Q&2 


bo  — 


h  = 


(ro  +  £*)3  £* 

fc4(3rg  +  3  r2£*  +  r0£*2)Q,2  2Qh 

(r0  +  £*)3  £* 

(3rg  +  3rgr  +  r0£*2)Q2 


+  4ff2, 


(ro  +  i*)3 


, .  (3rg  +  3 r2£*  +  3r0r2  +  £ *3  )02  , 
(h - 


(23) 

(24) 

(25) 


Sketch  of  proof:  The  linearization  of  (22)  upon  application  of  the  linear  state  feedback  T 
above  is  given  by 


where 


T+X  =  L°X 

dt 


0 

1 

0 

0 

0 

0 

\ 

-°! 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-a\ 

0 

0 

2fi 

e* 

0 

0 

0 

0 

0 

1 

0 

0 

-h 

-k2  +  2£*tt 

-k3  +  a3 

—  Aq 

/ 

Hence  the  characteristic  equation  of  the  closed-loop  system  is 


(A2  +  a2)  (A4  +  k\}?  +  6iA2  -(-  b2\  +  63)  —  0 


(26) 
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where  bi,  i  —  1,  3  are  as  defined  in  Eqs.  (23)-(25).  The  lemma  follows  readily  by 

applying  the  Routh-Hurwitz  test  to  the  second  polynomial  factor  in  the  left  side  of  Eq. 
(26). 

The  next  result  follows  readily  from  Lemma  1,  and  demonstrates  that  the  set  of 
feedbacks  of  Lemma  1  is  not  vacuous,  at  least  in  the  practically  interesting  case  r0  £*. 
The  result  also  holds  for  larger  £*,  but  the  corresponding  conditions  on  the  feedback  gains 
become  very  complicated. 

Corollary  1.  If  ro  £*,  the  conclusion  of  Lemma  1  holds  when  any  of  the  following  three 
conditions  is  satisfied: 

(i)  0  >  k\  >  l*Clk±(l  —  2ff0)  ^2  =  0,  /c3  >  3fl2,  and  k±  >  0; 

(ii)  0  >  kx  >  &2A4,  k 2  <  0,  >  3Q2,  and  &4  >  0; 

(iii)  k\  =  0,  &2  <  min{2flf*,  £3  +  3fl2},  k$  >  3f22,  and  >  0. 

From  the  closed-loop  characteristic  equation  (26),  we  infer  that  the  system  has  an 
uncontrollable  pair  of  pure  imaginary  eigenvalues.  Moreover,  it  is  easy  to  see  that  even 
when  the  states  (j>  and  are  used  in  the  tension  control  law,  these  two  pure  imaginary 
eigenvalues  remain  fixed.  Hence  the  stability  of  the  closed-loop  system  cannot  be  identified 
from  the  linearized  model  alone. 

The  closed-loop  system  (upon  application  of  a  feedback  law  as  in  Lemma  1  or  Corollary 
1)  is  approximated,  to  third  order  in  X,  by 

jtX  =  L0X  +  Qo(X,X)  +  C0(X,X,X)  (27) 

in  which  Lq  has  a  complex  conjugate  pair  of  pure  imaginary  eigenvalues,  with  the  remaining 
eigenvalues  in  the  left  half  of  the  complex  plane.  This  situation  is  an  example  of  a  critical 
case  in  nonlinear  stability  analysis.  Its  resolution  may  be  approached  via  results  on  Hopf 
bifurcation  for  one-parameter  families  of  nonlinear  systems  (see,  e.g.,  [16]).  The  local 
asymptotic  stability  of  the  origin  of  Eq.  (27)  can  be  concluded  from  the  negativity  of 
an  associated  “stability  coefficient,”  often  denoted  by  j32.  The  value  of  this  coefficient 
can  be  obtained  in  several  ways.  One  can,  for  instance,  study  normal  forms  of  Eq.  (27). 
Alternatively,  one  observes  that,  in  the  situation  at-hand,  smooth  parametrizations  of 
(27)  will  typically  exhibit  a  Hopf  bifurcation.  The  stability  of  the  bifurcated  periodic 
solutions,  as  well  as  that  of  the  origin  of  (27),  follows  from  the  negativity  of  the  Floquet 
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exponents  of  these  periodic  solutions.  The  stability  coefficient  /?2  can  be  obtained  as 
the  leading  coefficient  in  an  asymptotic  expansion  of  the  critical  Floquet  exponent.  The 
coefficient  /?2  may  be  computed  systematically.  When  /?2  <  0,  the  equilibrium  point  is 
locally  asymptotically  stable,  while  /?2  >  0  implies  instability  of  the  equilibrium.  The  case 
I32  =  0  is  inconclusive  regarding  stability.  We  now  proceed  to  use  an  algorithm  for  the 
computation  of  /?2  (see  for  instance  [10],  [16])  to  determine  the  dependence  of  /?2  on  the 
gains  ki,  i  =  1, . . .  ,4. 

Denote  by  r  and  l  the  right  (column)  and  left  (row)  eigenvectors,  respectively,  of  L0 
corresponding  to  the  imaginary  eigenvalue  ia\ .  Requiring  Ir  =  1,  we  have 

r  =  (1 ,  icq ,  0, 0, 0,  Of  (28) 

/  =  (^"i2r,0,0,0’°)  (29) 


Next,  solve  the  equations 


-L0a  =  -Qo(r,r), 


(2ia\I  —  Lo)b  =  -Qo(r,  r) 


(30) 

(31) 


for  the  vectors  a  and  b.  Here,  an  overbar  denotes  complex  conjugation  and  I  denotes  the 
identity  matrix.  We  find 

(32> 

b  —  (0,0,ci,c2,c3,c4)t  (33) 


where  the  a*  are  as  in  Appendix  A,  and  the  ct  and  dt  are  given  by 


.  a5  i*d\ 


ci  = 


d\  zc?2  2 

c2  =  2ia\Ci 


+  V-  --r(*3-a3)--w*} 


t 

C3  =  T- 


i£*(4af  —  a|) 


ci 


c  4 


4aiD 

it*a\  ^*(4af  —  a|) 


+ 


2Q 


ci 


(34) 

(35) 

(36) 

(37) 
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and 


di  =  h  +  ^(4«i  -  4) 

d2  =  2 ai(k2  -  2 £*S2)  +  r(4p  ^\a3  +  4a?  -  k3 ) 

*±(i\  S Z 


(38) 

(39) 


In  general,  the  value  of  the  bifurcation  stability  coefficient  j32  is  given  by  the  formula 


/32  =  2Re{2lQ0(r,  a)  +  lQ0(r,b )  +  -lC0(r,r,r)} 


(40) 


In  our  case,  we  then  have 


(6r%  +  4r%£*  +  r0£*2)Q2 


h  =  -{-ftlm(c2)  +  2(n+e.y 


a  i 


Jm(c3)  +  —  £e(c4)} 


4 

ai 


-Re(ci) 


(41) 


where  (assuming  £*  <  r0) 


j  o  ^*(4a? -a?)(6f3 +4r2r  +  r0r2)^2  ax(4a? -al) 

d3--2aiSl  (r0+^)4  + - 2ft - >0  (42) 


*(cO  =  d2  ^((y  +  ^  _  j(*»  -  «»))<*!  -  |tlOlClfc} 


(43) 


Since  a\  >  0,  d3  >  0,  and  since  stability  is  implied  by  ft2  <  0,  we  have  the  following 
result. 

Theorem  I.  If  a  linear  state  feedback  controller  T  as  defined  in  Lemma  1  is  applied,  with 


P*n 2  P*  1 

(f  +  -  j(*3  -  a3))di  -  ~Uaxtd2  <  0, 

then  the  equilibrium  point  xq  will  be  rendered  asymptotically  stable  for  the  system  (10)- 

(15). 

The  next  result  readily  follows  from  Theorem  1  and  Corollary  1. 

Corollary  2.  If  ?'o  ^  £*  and  a  linear  state  feedback  control  T  as  in  Lemma  1  is  applied, 
with  either  of  the  following  three  conditions  satisfied,  then  the  conclusion  of  Theorem  1 
holds: 
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(i)  0  >  k\  >  £*£tkffl  —  2!^)?  ^2  =  0,  3fi2  <  <  14a2,  and.  k 4  >  0; 

(ii)  0  >  Jbi  >  k2h,  0  >k2>  2£*n  -  i^(19  -  fa),  302  <  <  Hfi2,  and  h  >  0; 

(iii)  h  =  0,  min  {20T,3a2}  >  k2  >  2ttt  -  ±^(19  -  fa),  3 O2  <  <  14a2,  and 

&4  >  0. 

Case  2;  Nonlinear  State  Feedback 

Next  we  present  a  result  on  stabilization  with  a  tension  control  law  including  both 
linear  and  nonlinear  terms.  The  nonlinear  terms  introduce  more  flexibility  in  the  design, 
and,  as  will  be  seen  in  Section  IV,  can  lead  to  superior  transient  response. 

By  computing  the  eigenvectors  l  and  r  and  using  the  formula  (40)  for  ft2 ,  we  find  that 
any  cubic  term  in  the  feedback  has  no  effect  on  the  value  of  ft2.  We  are  therefore  led  to 
hypothesize  a  feedback  containing  only  linear  and  quadratic  terms.  The  component  of  the 
closed-loop  quadratic  term  Q o(X,X)  depending  on  the  states  6,  toe,  t  or  v  is  also  found  to 
have  no  effect  on  ft 2.  The  next  theorem  gives  conditions  for  a  nonlinear  feedback,  of  the 
form  motivated  by  these  observations,  to  be  stabilizing. 

Theorem  2.  If  condition  (44)  below  holds,  then  the  applied  tension  control  force  T  = 
m(—U  —  ki6  —  k2u>o  —  k^£  —  k±v  —  q\ ft2  —  q2<f>to<f>  —  stabilizes  the  system  (10)-(15), 
where  the  k i ,  i  =  1 ,  . . . ,  4  satisfy  the  conditions  of  Lemma  1 . 

(  gl2~a--  +  ~  “  a3)  +  y(“92  -  £*h)d2  <  0  (44) 


Here,  d\,  d2  are  given  in  Eqs.  (38)  and  (39),  and  the  a*  are  as  specified  in  Appendix  A. 


The  proof  entails  checking  the  effect  on  the  value  of  ft2  of  adding  the  extra  quadratic 
term  -(qift2  +  q2ftw<f)  +  in  the  last  row  °f  Q o(X,X). 

As  mentioned  above,  inclusion  of  nonlinear  terms  in  the  feedback  control  may  be  used 
to  improve  the  transient  response  of  the  stabilized  system.  In  particular,  the  rate  at  which 
system  trajectories  decay  toward  the  equilibrium  point  may  be  significantly  increased. 
Simulation  evidence  for  this  is  given  in  Section  IV. 

It  is  not  difficult  to  give  analytical  reasoning  to  support  this  conclusion,  and  to  guide 
in  the  tuning  of  the  linear  and  quadratic  feedback  gains.  Assume  r0  >  £*,  and  use  Eq. 
(40)  with  a  feedback  of  the  form  given  in  Theorem  2  to  ascertain  the  approximate  formula 


g2  „ - - - {(-^ 

12  ~  2 (d2  +  d2y{  2 


1 


+  - 


£* 

(h  ~  a3)—)d 1 
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(45) 


+  ~^(—Q2  —  £*ki)d2} 

Eq.  (45)  can  be  used  to  show  that,  if  r0  £*  and  the  linear  gains  ki ,  *  =  4,  are 

chosen  according  to  condition  (i)  of  Corollary  2,  with  k\  =  0,  then  j32  may  be  rendered  as 
negative  as  desired  simply  by  setting  the  quadratic  gains  q2  =  q%  =  0  and  taking  q\  >  0 
and  sufficiently  large.  Thus  the  gains  ki  may  be  used  to  place  four  of  the  eigenvalues  of 
(10)-(15)  in  the  left  half  of  the  complex  plane,  while,  independently ,  the  gains  qi,  i  =  1,2,3, 
are  used  to  make  j32  negative  and  of  large  magnitude. 

III.3.  Stabilization  for  In-Plane  Angle  Near  9  =  tt 

Similarly,  now  let  x„  =  (0, 0,  tt,  0,  £*,0)T  denote  the  equilibrium  point  of  interest,  and 
X  =  x  —  Xn  be  the  differential  state  variation.  Then  the  system  (10)-(15),  to  third  order 
near  the  equilibrium  point  x„,  may  be  written  as  follows 


d_ 

dt 


X  =  L„X  +  Q„(X,  X)  +  Cx(X,  X,  X)  +  eUn  +  e 


T 

m 


Here,  Lx,  Q C ^  are  as  identified  in  Appendix  B.  The  next  lemma  is  analogous  to  Lemma 
1,  and  so  is  stated  without  proof. 

Lemma  2.  Let  the  applied  tension  force  be  of  the  form  T  =  m(— Un  —  k\9  —  k2L0g  —  k^i— k<iv), 
where  6  :=  6  —  n  and  u>g  6  =  9.  Then  the  “in-plane”  Jacobian  matrix  of  Eqs.  (12)-(15), 

i.e.,  the  Jacobian  of  the  right  side  of  (12)-(15)  with  respect  to  (9,ojg,£,  v ),  will  be  stable  at 
the  equilibrium  point  xn,  if  ki,  i  =  1,  . . .,  4,  satisfy  the  following  conditions: 


(i)  ki  >  0, 

(ii)  h\,h2,hz  >  0,  and 


(iii)  kihih2  —  h\  —  k2h$  >  0. 

Here,  the  auxiliary  parameters  hi,  h2,  h3  are  given  by 


hi  =  k3  — 


(2r0l*2  —  l*3)fl2  _  2 tlk2 
(r0  -  £*)3  £* 


+  402, 


jfc4(3rg  -3r2£*  +  r0f2)fi2  2 
h2  ~  (r0  -  £*f  L*  ’ 

,  (3 r3  -  3r2£*  +  r0£*2)£l2  n  (3rg  -  3r2^  +  3r0£*2  -  £*3)Q,2  N 

h3~  (r0  —  £*)3  (3  (r0  —  £*)3  } 


(46) 

(47) 

(48) 
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Corollaries  1  and  2  remain  valid.  The  detailed  statements  need  not  be  given. 

IV.  Simulation  Results 

A  TSS  with  the  following  characteristics  is  considered: 

•  Nominal  tether  length  l*  —  100  km, 

•  Orbital  radius  ro  =  6598  km, 

•  Subsatellite  mass  m  —  170  kg, 

•  Orbital  angular  velocity  —  0.0011781  radians /second. 

Let  the  equilibrium  point  of  interest  of  (10)-(15)  be  xq  —  (0,0,0,0,£*,0).  Simulation 
results  will  now  be  presented  which  illustrate  the  system  dynamics  for  the  various  types 
of  control  studied  in  this  paper. 

Let  the  initial  conditions  of  the  system  be  (f>  =  0.01  radians,  0  =  —0.01  radians,  and 

toe  =  Uj,  =  0. 

Example  IV.  1.  (  No  tension  control:  latch  mechanism) 

Suppose  the  reel  mechanism  acts  like  a  latch  fixing  I  at  I*.  The  system  response 
for  the  in-plane  angle  0  and  the  out-of- plane  angle  <j>  is  shown  in  Figures  3(a)  and  3(b), 
respectively.  We  observe  an  apparent  undamped  oscillation  near  the  equilibrium  point  xo . 

Example  IV.  2.  (Linear  stabilizing  feedback) 

The  tension  controller  is  taken  as  T  =  —m{U  +  k$I  +  k±v),  with  k$  —  3.1fl2,  — 

0.0034,  and  U  —  0.41019.  The  control  law  is  stabilizing,  as  can  be  checked  using  Theorem 
1.  Indeed,  02  —  —0.0004  for  the  closed-loop  system.  The  response  of  the  variables  of 
0,  and  the  deviation  I  of  the  tether  length  are  shown  in  Figures  4(a),  4(b)  and  4(c), 
respectively.  However,  it  is  not  easy  to  see  in  Figure  4(a)  any  decay  of  the  oscillation  in 
the  out-of-plane  angle  <f>.  This  may  be  attributed  to  the  fact  that  \02  \  is  small.  The  applied 
tension  force  is  shown  in  Figure  4(d). 

Example  IV. 3.  (Linear-plus- quadratic  stabilizing  feedback) 

Let  the  tension  control  law  be  of  the  form 

T  =  —m(U  +  hi  +  k±v  +  qi<j)2  +  q2<t»x>4>  +  qz<^>%),  (49) 
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where  U  =  0.41019.  The  out-of-plane  angle  <j>  decays  when  fc3  =  3.1172,  Aq  =  0.0034, 
qi  =  1500,  and  q2  =  qz  —  0  as,  shown  in  Figure  5(a).  However,  this  is  at  the  expense  of 
large  variations  in  9  and  £,  as  depicted  in  Figures  5(b)  and  5(c).  The  applied  tension  force 
is  shown  in  Figure  5(d). 

Example  IV. 4-  (Linear-plus- quadratic  stabilizing  feedback) 

A  further  example  is  depicted  in  Figure  6.  In  this  example,  &:3,  k4  and  <?3  are  unchanged 
from  their  previous  values  (given  in  Example  IV. 3),  but  now  qi  =  0,  and  q2  =  106. 

Example  IV. 5.  (Switching-type  stabilizing  feedback) 

Figure  7  relates  to  an  example  invoking  a  switching  control  strategy.  The  nonlinear 
feedback  control  law  of  Example  IV. 3  is  used  for  the  first  5  hours  of  the  simulation.  Then 
the  control  law  is  switched  to  a  purely  linear  feedback  with  the  parameters  values  specified 
in  Example  IV. 2. 

V.  Concluding  Remarks 

In  this  paper,  we  have  presented  analytical  designs  of  tension  feedback  control  laws 
for  the  stabilization  of  the  tethered  satellite  system  during  station-keeping.  These  designs 
are  based  on  calculations  related  to  Hopf  bifurcation  stability.  The  calculations  have  been 
performed  for  a  model  of  the  tethered  satellite  system  derived  under  several  simplifying 
assumptions.  This  model  is  characterized  by  its  nonlinearity  and  the  existence  of  two 
critical  modes,  one  of  which  cannot  be  removed  by  (linear)  feedback.  Notwithstanding  this 
fact,  we  have  been  able  to  construct  stabilizing  controllers  using  linear  and/or  quadratic 
feedback.  Cubic  terms  were  not  included  in  the  feedback  laws  since  the  nonlinear  stability 
calculations  indicated  that  their  effect  might  be  of  only  secondary  significance.  Moreover, 
simulation  was  used  to  demonstrate  the  validity  of  the  analytical  designs.  The  simulations 
also  indicated  the  importance  of  quadratic  feedback  of  the  out-of- plane  angle  <j>  and/or 
the  out-of-plane  angular  rate  in  improving  the  transient  response  of  the  out-of-plane 
variables,  i.e.,  in  dampening  the  roll  oscillations. 

Regarding  the  issue  of  how  to  achieve  further  improvements  in  the  transient  response, 
several  possibilities  arise.  Optimization-based  computer-aided  design  tools  can  be  applied 
to  systematically  search  for  linear  and/or  nonlinear  control  gains  resulting  in  a  subopti- 
mal  transient  response.  If  other  actuators,  such  as  subsatellite  thrusters  or  tether  base 
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movement,  are  available  in  addition  to  tether  tension  control,  then  one  expects  improved 
transient  performance. 
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Appendix  A 

The  values  of  the  coefficients  a,-,  i  =  1,  . . .,  15  are  as  listed  below. 
,(4r0J+6r02r  +  4r„C2+C3)u/2„ 

“1=( - WTW - >  n 


a  2 


«3 


_  ,(3rg+3rgr+r0r2)u/2o 
(  (ro+£*)3  ’ 

(3  rg  +  3  r2l*  +  3r0l*2  +  l*3)f l2 
(r0  +  i*)3 


(6rgf  +  6rj)£*2  +  4  r2£*3  +  r0£*4)Q2 
Q4~~  2  (ro-M*)4 

(8  r$£*  +  14r^*2  +  1 6rjj£*3  +  9r0^*4  +  2^*5)02 
°5  2 (r0+£*)4 
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(6rp  +  9  r%£*  +  lOrgl*2  +  5r2l*3  +  r0£*4)Q2 
2  (r0  +  £*)5 

(16rg  +  29r$£*  +  50  r$£*2  +  45r2l*3  +  21  r0£*4  +  4  £*5)fl2 

6(ro  +  £*)5 

(10r3  +  5  r20£*  +  r0£*2)Q2 
(r0  +  £*)5 

(12  r3  +  9  r4£*  +  10r3l*2  +  5  r2l*3  +  r0l*4)O2 

6(ro  +  £*)5 


«io 


an 


(27r4l*  +  30r[jl*2  +  15  rgl*3  +  3r0l*4)f22 
6(r0  +  £*)5 


(4rp  +  2  r\t  +  10  r3l*2  +  lOrgl*3  +  5  r0£*4  +  £*5)Q2 
(r0  +  £*)5 


a\2 


(6r|  +  4r2l*  +  r0£*2)Sl2 
(r0  +  ^*)4 


3r3ft2 

ai3_  (r0  +  £*)4 


ai4  = 


(3rp  -  3r4l*)02 
(r0  +£*)5 


_  4r3H2 

(r0  +  £ *)5  ' 


Appendix  B 

The  system  model  (10)-(15)  is  approximated,  to  third  order  in  the  states,  near  the 
equilibrium  point  xn  by 


X  =  LnX  +  QW(X,  X)  +  C^X,  X,  X)  +  eUn  +  e— 


where 

1  0  0  0  0  \ 

0  0  0  0  0  ’ 

0  0  10  0 
0  -/|  0  0  -f- 

0  0  0  0  1 

0  0  2ffi  /3  0  / 


/  0 

0 
0 
0 

\  0 
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-20(^o >o  +  -  j?LO<j>v 

0 

f\2@£  —  JZ’UqV  +  2Q,(f)U><j>  +  j^-iv 

0 

\fj2  +  £*u2  +  2W  +  /5</*2  +  £*u2  +  hie  / 


f§02<t>  —  +  /V^3  +  /8^2  +  ~jfc2U<l>$V 

0 

f903  +  /io<¥2  +  /s^2  +  2  fa+wo  +  -jkueiv  -  $£2v 

0 

V  -  2t£luje<$>2  +  hift  +  +  /15^3  / 


e  =  (0,0, 0,0, 0,1)' 

_  (3rgl*  -  3rpl*2  +  l*3)02 
(r0  -  f?*)2 


and  the  values  of  /;,  i  =  1,  . . .,  15  are 


„  A4rZ-6r2t*+4ror2-£*3),1/2„ 

/l=(  (r0-£*)3  j 

,  A3r30-Zr2£*+r0£*2),1/2„ 

h  =  {  (^F?  } 


h 


h  = 


/«  = 


U 


h 


h  = 


(3r3  -  3r2l*  +  3r0l*2  -  l*3)ft2 
(r0  -^*)3 

(6rj*l*  -  6r3l*2  +  4r2l*3  -  r0£*4)f22 
2(r0  -  £*Y 

(8r|T  -  14r3l*2  +  16r2l*3  -  9r0l*4  +  2l*5)^2 
2(r0  —  £*)4 

(6rg  -  9 r%£*  +  10r(|f  2  -  5 r%£*3  +  r0l*4)Q2 
2(r0-^)5 

(16r|  -  29r4l*  +  50r3l*2  -  45r2l*3  +  21r0l*4  -  4l*5)ft2 

6(r0  -  £*f 

(10r3  -  5 r2£*  +  r0£*2)Sl2 
(r0  -  £*f 
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fa  = 

/io  = 
/ll  = 
/l2  = 
/l3  = 
/l4  = 
/l5  = 


(12r|  -  9 rgl*  +  lOrgl*2  -  5r2l*3  +  r0l*4)ft2 
6(ro  —  ^*)5 

(27rgl*  -  30rj]l*2  +  15r2l*3  -  3 r0l*4)ft2 
6(r0  -^*)5 

(4rp  -  2rjjl*  +  10 r%£*2  -  IQrgl*3  +  5 r0l*4  -  l*5)ft2 

(r0  -£*)5 

(6r3  -  4 rgl*  +  r0l*2)Q2 
(r0  -  ^*)4 

3rpfl2 

(ro-^)4 

(3rp  +  3rpl*)n2 

(r0  -^*)5 

4rpf22 
(r0  -£*)5' 


Figure  1.  TSS  in  circular  orbit. 
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1/ 

—(c)  t  (hr.)  0  •  ( d )  4  •  t  (hr.) 

Figure  5.  Simulation  results  for  nonlinear  feedback  system  ( q\  =  1500). 
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